Entanglement Entropy as a Portal to the Physics of Quantum Spin Liquids 
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Quantum Spin Liquids (QSLs) are phases of interacting spins that do not order even at the 
absolute zero temperature, making it impossible to characterize them by a local order parameter. 
In this article, we review the unique view provided by the quantum entanglement on QSLs. We 
illustrate the crucial role of Topological Entanglement Entropy in diagnosing the non-local order in 
QSLs, using specific examples such as the Chiral Spin Liquid. We also demonstrate the detection 
of anyonic quasiparticles and their braiding statistics using quantum entanglement. In the context 
of gapless QSLs, we discuss the detection of emergent fermionic spinous in a bosonic wavefunction, 
by studying the size dependence of entanglement entropy. 
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I. INTRODUCTION 

Quantum Spin Liquids (QSLs) are phases of strongly 
interacting spins that do not order even at absolute 
zero temperature and therefore, are not characterized 
by a Landau order parameter—. They are associated 
with remarkable phenomena such as fractional quantum 
numbers^Ti^, transmutation of statistics (eg. fermions 
appearing in a purely bosonic model)^ii, and enabling 
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Figure 1: Following Ref. 1191 , consider four regions A1-A4 
(shaded cyan). The thickness as well as the linear size of 
all regions Ai is much bigger than the correlation length ^ of 
the system. The combination {Si — S2) — (S3 — 5*4) equals 
zero for short-range entangled systems. To see this, we note 
that Si — S2 and S3 — S4 correspond to the change in en- 
tanglement entropy associated with closing the region A2 and 
A4, respectively, at the top. For a wavefunction with short- 
range entanglement, the difference (Si — S2) — (S3 — S4) should 
vanish as e~^^^ where R is approximately the linear size of 
regions Ai. Therefore, a non-zero 7 implies that the system 
is long-range entangled. 



otherwise impossible quantum phase transitions^, to 
name a few. Spin liquids may be gapless or gapped. 
While current experimental candidates for spin liquids 
appear to have gapless excitations^, gapped spin liquids 
are indicated in numerical studies on the KagomeiSiii, 
squarei^ii^ and honeycomb lattice^. A distinctive prop- 
erty of gapped spin liquids is that they are characterized 
by topological order - i.e. ground state degeneracy that 
depends on the topology of the underlying spac o^^'^^ . 
The gapless spin liquids are relatively difficult to char- 
acterize and their low-energy description often involves 
strongly interacting matter-gauge theoriesii. 

In the theoretical search for QSLs, two qualitatively 
different kinds of questions can be asked: 
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• Given a particular Hamiltonian H, does it admit 
a QSL as its ground-state? If yes, can one further 
characterize the QSL? 

• Given a ground state wavefunction(s), can one tell 
if it corresponds to a QSL? If yes, can one further 
characterize the QSL? 

These two questions are related: if a numerical scheme 
can provide one with access to approximate ground state 
wavefunction of a particular Hamiltonian iJ, then an- 
swering the second question leads to an answer to the 
first one as well. The recent years have seen remark- 
able progress in tackling these questions. Of particular 
usefulness has been the application of ideas from quan- 
tum information science to quantum many-body physics. 
Motivated by these developments, in this review we will 
focus on the second question above from an information 
theoretic point of view, and we will discuss some of the 
advances on the first question in the section IVIII 

The usefulness of entanglement entropy to understand 
the physics of QSLs becomes apparent when one intro- 
duces the notion of 'short-range' and 'long-range' entan- 
glement. The heuristic picture of a short-range entangled 
phase is that its caricature ground state wavefunction can 
be written as a direct product state in the real space, 
hence the nomenclature. The phases of matter that do 
not satisfy this property are thereby called 'long-range 
entangled'. From this definition, gapless phases, such as 
a Fermi liquid or superfluid, that exhibit algebraic de- 
cay of local operators are always long-range entangled in 
the real space. Naively, one might think that the con- 
verse would hold, that is, existence of short-range corre- 
lations for local operators will imply short-range corre- 
lation. It is precisely here that the connection between 
gapped QSLs and ground state entanglement properties 
comes into play. To wit, gapped QSLs are phases of mat- 
ter that show long-range entanglement in their ground 
states despite the fact that the correlation functions of 
all local operators decay exponentially with distance. 

To make the connection between gapped QSLs and 
long-range entanglement precise, let us define some ba- 
sic terminology that will also be useful throughout this 
review. Given a normalized wavefunction, |(/)), and a par- 
tition of the system into subsystems A and B, one can 
trace out the subsystem B to give a density matrix on A: 
Pa = TrB|(/))(0|. The Renyi entropies are defined by: 

Sn = ^\0g{TTpl) (1) 

1 — n 

It is common to pay special attention to the von Neu- 
mann entropy, 5*1 = —l^v[p A^og p a] (obtained by taking 
the limit — > 1 in Eqn[T]). However, the Renyi entropies 
are also equally good measures of entanglement, and are 
often easier to compute. 

In seminal papers by Kitaev, Preskilli^ and Levin, 
Werti^, it was shown that the entanglement entropy of 
a topologically ordered phase, such as a gapped QSL, as 



a function of the linear size of disc-shaped region A of 
linear size Z, receives a universal subleading contribution: 

^(L) = --6o7 + 0(l//) (2) 
a 

where 6o is the number of connected components of the 
boundary of region A and the constant 7 is a univer- 
sal property of the phase of matter that is non-zero if 
and only if the phase is topologically ordered. One may 
wonder why does a non-zero 7 implies long-range entan- 
glement? One way to understand thisi^ is to consider 
the combination Sq = {Si — S2) — {S3 — Si) for the geom- 
etry shown in FiglT] If the phase was in fact short-range 
entangled, then Sq will vanish identically. From Eqn|5J 
So = —2"f and therefore, a non-zero 7 implies long-range 
entanglement. A different way to derive the same conclu- 
sion is to define a quantity Siocai such that it is obtained 
by locally patching the contributions to the entanglement 
entropy from the boundary OA of region A^. If the phase 
is short-range entangled, then Siocai is the entire S. As 
shown in Ref. [IS, on general grounds, Siocai cannot con- 
tain a constant term, and therefore, 7 7^ implies long- 
range entanglement. 

Gapless spin-liquids pose a different set of questions. 
The gapless excitations make it difficult to assign them 
a simple topological number such as 7 above. One ex- 
ception is gapless spin-liquids where the emergent gauge 
degrees of freedom are gapped and are coupled to gap- 
less matter. In this case, there are known contributions 
to the entanglement entropy that can characterize the 
nature of QS L^^'^^ . In this review, we will focus on cases 
where the emergent gauge fields are instead gapless. As 
an example, consider a two-dimensional QSL where the 
gapless matter fields are neutral fermions which form a 
Fermi surface, and are coupled to a gapless gauge field. 
As we discuss in Sec lVIi empirically one finds that their 
entanglement entropy scales as S" ~ Hog I, akin to gap- 
less fermions in two dimensions^i^, thus suggesting the 
presence of emergent fermions in a purely bosonic state. 

II. ENTANGLEMENT ENTROPY AND MANY 
BODY PHYSICS 

One of the most interesting properties of EE is that 
S{A) = S{A) where A is the complement of region A. 
This suggests that EE is a property of the boundary dA 
of region A since dA = dA, which further suggests that 
for a short-ranged correlated system, EE scales linearly 
with the the boundary dA: S{1) ^ /''"^ where / is the 
linear size of region A^^^. Such a scaling of EE goes by 
the nomenclature of "Area Law" or "Boundary Law" . In- 
deed, for one-dimensional systems with a gap, this state- 
ment can be proved rigorously^ and it has also been ob- 
served numerically for various gapped two-dimensional 
systems^^. Furthermore, even when the correlations of 
local operators arc long-ranged, for example, in a free 
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fermionic system with a Fermi surface, the Boundary 
Law is violated by only by a multiplicative logarithm, 
that is, S{1) - l'^-^ logi^l'p.^, while for a conformal field 
theory (CFT), where the low-lying modes exist at only 
discrete point(s) in the momentum space, the Boundary 
Law is not violated at all in d > 1—. 

To give a flavor of the information theoretic approach 
to many-body physics, let us consider one of the earli- 
est such applications that concerns the relation between 
the central charge c of a one-dimensional CFT, and the 
coefficient of the logarithmic term in the bipartite entan- 
glement entropy S for the corresponding CFT—. The 
result is 

5(/) = ^log(-) (3) 
D a 

where I is the size of the subsystem A and a is a mi- 
croscopic cutoff. The main observation to gather from 
this equation is that just given the ground state wave- 
function of a CFT, its central charge can be extracted by 
calculating the entanglement entropy. Note that rescal- 
ing where a is a non- universal number does not 

change the coefficient of the logarithm. This makes phys- 
ical sense, since c is a property of the low-energy theory 
and should be insensitive to the details of how the theory 
is regularized at the level of the lattice. 

An extension of this result to a two dimensional CFT 
isSS: 



s{i) = --r (4) 

a 

for a region A with linear size I and smooth boundaries. 
Here F is a constant that depends only on the shape of 
region A. Rescaling I ^ l/a does not change F akin to 
the case of one-dimensional CFT, albeit in this case F 
will generically be shape dependent. As we will see, this 
result will be useful to characterize gapless QSLs whose 
low-energy theory is a CFT (Sec lVI|) . 

The common feature of the above two examples is that 
EE as a function of system size / contains terms that 
are universal properties of the low-energy description of 
the system. Returning to our topic of interest, namely 
QSLs, one of the most interesting resul t "'^^'^^ is that for 
gapped systems in two dimensions, the subleading term 
in EE is universal and contains information regarding 
the presence or absence of topological order, and also 
partially characterizes the topological order: 

^(0 = --7 + 0(1/0 (5) 
a 

The subleading constant 7 is called "Topological En- 
tanglement Entropy" (TEE) and as already mentioned in 
the introduction, is a universal property of the phase of 
the matter, even though one happens to have has access 



to only one of the ground states in the whole phase dia- 
gram. 7 is non-zero if and only if the system has topolog- 
ical order. Furthermore, 7 = log(D) where D = y^iW 
is the so-called total quantum dimension associated with 
the phase of matter whose quasiparticles carry a label 
di, the individual quantum dimension for the z'th quasi- 
particle. As a concrete example, for Kitaev's Toric code 
model^, there are four quasiparticles, the identity parti- 
cle 1, the electric charge e, the magnetic charge m, and 
the dyon em. The quantum dimension is one for each 
of these quasiparticles and hence the total quantum di- 
mension D = \l\ = 2. We note that the the universal 
contribution to EE was first observed in the context of 
discrete gauge theories in Ref. 

The Eqn[S] serves as a definitive test to see if a par- 
ticular ground state corresponds to topological ordered 
state or not. If the answer is indeed in affirmative and 
if the state also possesses a global SU{2) spin-rotation 
symmetry, then it is a gapped QSL. Having established 
that the wavefunction corresponds to QSL, it is imper- 
ative to further characterize the state further since the 
total quantum dimension D only a partial characteriza- 
tion of topological order. For example, for D = 2, there 
exist two distinct QSLs, the Z2 QSL, which has the same 
topological order as the Kitaev's Toric code^, and the 
doubled semionic theory^, which can be thought of as a 
composite of Laughlin v = 1/2 quantum Hall state and 
its time-reversed partner. To achieve a detailed charac- 
terization of topological order, TEE again turns out to be 
a very useful concept. As we will show in Sec lVB ll TEE 
can be used to extract not only the individual quantum 
dimensions di for the quasiparticles, but there braiding 
statistics as well. 

An elegant different approach to a more complete iden- 
tification of topological order is through the study of the 
entanglement spectrum^. This method is extremely use- 
ful for phases that have a gapless boundary such as quan- 
tum Hall and topological insulators/superconductors, 
though it may not be applicable for topological phases 
such as the Z2 spin liquid that do not have edge states. 



III. VARIATIONAL WAVEFUNCTIONS FOR 
QUANTUM SPIN LIQUIDS 

A quantum spin liquid (QSL) can be viewed as a state 
where each spin forms a singlet with a near neighbor, but 
the arrangement of singlets fluctuates quantum mechan- 
ically so it is a liquid of singletsii^. Theoretical models 
of this singlet liquid fall roughly into two categories. In 
the first, the singlets are represented as microscopic vari- 
ables as in quantum dimer and related models^"—, and 
are suggested by large N calculations^. Topological or- 
der can then be established by a variety of techniques 
including exact solutions. In contrast there has been less 
progress establishing topological order in the second cat- 
egory, which are SU(2) symmetric spin systems where 
valence bonds are emergent degrees of freedom. More 
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than two decades ago, Andcrson^i proposed constructing 
an SU(2) symmetric spin liquid wavefunction by starting 
with a BCS state, derived from the mean field Hamilto- 
nian: 

H=- Y,{trr'flrU,r' + Arr' fljl,,} + h.C. (6) 
rr' 

and Gutzwiller projecting^ it so that there is exactly one 
fermion per site, hence a spin wavefunction. Variants of 
these are known to be good variational ground states 
for local spin Hamiltonians (see e.g4ii^) and serve as 
a compact description of most experimental and SU{2) 
symmetric liquids. Approximate analytical treatments 
of projection, that include small fluctuations about the 
above mean field state, indicate that at least two kinds of 
gapped spin liquids can arise: chiral spin liquids^ (CSL) 
and ^2^'^'^ spin liquids. The low energy theory of CSL 
corresponds to a Chern-Simons theory while in the Z2 
spin-liquid, it corresponds to dcconfined phase of the Z2 
gauge theory. However, given the drastic nature of pro- 
jection, it is unclear if the actual wavefunctions obtained 
from this procedure are in the same phase as the ground 
state of a Chern-Simons theory/Z2 gauge theory. In fact, 
detecting deconfinement in a gauge theory in the presence 
of matter-field is considered notoriously hard^. As we 
will see below, quantum entanglement provides an ele- 
gant solution to this problem. 

In contrast to gapped QSLs, the low-energy description 
of gapless QSLs consist of strongly interacting matter- 
gauge theories which are harder to analyze^i. Remark- 
ably, the above procedure of constructing mean-field 
spin-liquid states and using Gutzwiller projection to ob- 
tain wavefunctions can be applied to obtain gapless QSLs 
as well. To do so, one sets the BCS gap Arr' in EqnlB] 
to zero. This Gutzwiller projection technique is known 
to work well in one dimension where the projected Fermi 
sea spin wavefunction captures long distance properties 
of the Heisenberg chain, and is even known to be the 
exact ground state of the Haldane-Shastry^ Hamilto- 
nian. Though similar rigorous results are not available 
in two dimension, a critical spin liquid, the spinon Fermi 
sea (SFS) state, has been invoked^i^ to account for 
the intriguing phenomenology of aforementioned trian- 
gular lattice organic compounds^. In the SFS state, the 
spinous / hop on the triangular lattice sites giving rise 
to a Fermi sea, while strongly interacting with an emer- 
gent 'electromagnetic' t/(l) gauge field. The metal like 
specific heat and thermal conductivity seen in these ma- 
terials is potentially an indication of the spinon Fermi 
surface. In this regard, the Gutzwiller projected Fermi 
sea is known to have excellent variational energy for the 
J2 — Ji spin model on the triangular lattice, which is 
believed to be appropriate for the aforementioned trian- 
gular lattice compounds^. 

As mentioned in the introduction, the main task of 
this review is to illustrate the extraction of the universal 
properties of a QSL from the ground state wavefunction. 



Therefore, we now turn to the technical challenge of cal- 
culating EE of Gutzwiller projected wavefunctions that 
are candidates for QSL. 

IV. FROM QSL WAVEFUNCTIONS TO 
ENTANGLEMENT ENTROPY 

In this section, using a Monte Carlo scheme, we out- 
line how to obtain the entanglement entropy of a ground 
state function that has the form of a Gutzwiller pro- 
jected Slater determinant. As mentioned in the previous 
sections, such wavefunctions are a natural candidate for 
QSLs. We will focus on the second Renyi entropy S2 (see 
Eqn. [T]) since it suffices to establish topological order and 
is easiest to calculate within our scheme. 

Let (j)(a, b) be the wavefunction of interest, where a (b) 
be the configuration of subsystem A (B). Then the re- 
duced density matrix is pa{ci, a') = 0*(a, b)(j){a' , b)/N 
where the normalization N — X^a 6 ^)'?^(^' "^^^ 
second Renyi entropy ^2 is: 

= ]^ E r{a,b)cf>{a',b)4,*{a\b')c^{a,b') (7) 

aa'hb' 

One way to interpret this expression is to con- 
sider two copies of the system, in configuration 
\a,b)\a' ,b'). Then, the product wavefunction |<I>) = 
J2a b a' b' 4'{0'i b)(f>{a' , b')\a, b)\a', b') appears in the expres- 
sion above. In addition if we define the Swapj^ opera- 
tor, following ReS^, which swaps the configurations of 
the A subsystem in the two co^\es:Swapj^\a,b)\a' ,b') = 
\a' , b)\a, b') then Renyi entropy can be written simply as: 

' " — m — ^' 

The above expression suggests a Variational Monte Carlo 
algorithm for the Renyi Entropy. 

Defining configurations ai — a, 6, 012 — a', b' and /3i = 
a', b, 132 ~ a, 6', and 4>ai = 4>{0',b) etc., Eqn.® can be 
rewritten in the suggestive form: 

f{ai,a2) (9) 

0102 

Here the weights pa^ = \4'ai\'^ /J2a- l^^il"^ '^O'^" 
negative and normalized. The quantity / (01,012) = 
t'^^t'^^ when averaged over Markov chains generated 
with the probability distribution p(ai)p(a2) yields 
{Swap a)- The Gutzwiller projected wavefunctions to de- 
scribe the QSLs we consider in this review can be written 
as products of determinants, as explained in the previ- 
ous section. Such wavefunctions can be efficiently evalu- 
ated, which forms the basis for Variational Monte Carlo 
technique^. 



State 


Expected 7 


/ 

/calculated/ /expected 


Unprojected (^u = 1) 







-u.uuuozt u.uuoy 


Chiral SL La=3 


log 


V2 


0.99 ± 0.03 


Chiral SL La =4 


log 


V2 


0.99 ± 0.12 


Lattice u — 1/3 


log 


^/3 


1.07± 0.05 


Lattice u — 1/4 


log 


\/4 


1.06 ± 0.11 



Table I: Topological Entanglement Entropy 7caicuiatGd of spin- 
liquid states obtained using Monte Carlo technique discussed 
in the main text. 7expected is the exact value expected in the 
limit ^/i — >■ where ^ is the correlation length and I is the 
subsystem size. 



As shown in Ref. [ssUs^ the above algorithm, with 
slight modifications that we will discuss briefly in Scc lVIl 
can be used to calculate S2 within a few percent of error 
bar for subsystems with linear size La ^ 10, embedded 
in systems with linear size L <20. 



V. GAPPED QSLS AND ENTANGLEMENT 

Model systems such as Toric code^ and quantum 
dimer models^ provide a testing ground for the valid- 
ity of EqnlSl More realistic interacting quantum systems 
pose a greater challenge since the analytical methods are 
limited. Some of the earlier numerical work to extract 
TEE included exact diagonalization studies on small sys- 
tems, looking at quantum Hall Laughlin states^ and per- 
turbed Kitaev toric code models^. In the context of 
QSLs, a recent quantum Monte Carlo study^ on a sign 
problem free Hamiltonian used TEE to detect Z2 topo- 
logical order. This was a positive definite wave-function, 
with U(l) rather than SU(2) spin symmetry. In the rest 
of this section, we will focus on SU{2) symmetric spin- 
liquids^, using the technique explained in the previous 
section. 



A. Establishing Topological Order in Gapped QSLs 

In this subsection we discuss few examples where using 
the algorithm to calculate the entanglement entropy dis- 
cussed in the previous section, one can establish topologi- 
cal order in SU{2) symmetric QSLs^. As a first example, 
let us consider an SU (2) symmetric lattice wavefunction, 
the Kalmeyer-Laughlin Chiral Spin Liquid (CSL). 

The Chiral SL is a spin SU(2) singlet ground state, 
that breaks time reversal and parity symmetr y "^^i^^ . A 
wavefunction in this phase wave function may be ob- 
tained using the slave-particle formalism by Gutzwiller 
projecting a d + id BCS stated. Alternately, it can be 
obtained by Gutzwiller projection of a hopping model on 
the square lattice. This model has fermions hopping on 
the square lattice with a tt fiux through every plaquctte 
and imaginary hoppings across the square lattice diago- 
nals: 
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Figure 2: Illustration of a square lattice hopping model con- 
nected with a d+id superconductor. While the nearest neigh- 
bor hopping is along the square edges with amplitude t {—t 
for hopping along dashed lines), the second nearest neigh- 
bor hopping is along the square diagonal (arrows in bold), 
with amplitude -l-iA (— iA) when hopping direction is along 
(against) the arrow. The two sublattices in the unit cell are 
marked as A and B. 



H = Y}^,l}l,+iY,A,,fUk (10) 

Here i and j are nearest neighbors and the hopping 
amplitude is t along the y direction and alternating 
between t and —t in the x direction from row to row; and 
i and k are second nearest neighbors connected by hop- 
pings along the square lattice diagonals, with amplitude 
Aife = zA along the arrows and A.^fc — — iA against the 
arrows, see Fig. 1. The unit cell contains two sublat- 
tices A and B. This model leads to a gapped state at 
half filling and the resulting valence band has unit Chern 
number. This hopping model is equivalent to s. d + id 
BCS state by an SU (2) Gauge transformation^. We use 
periodic boundary conditions throughout this section. 

Due to the fact that this Hamiltonian contains only 
real bipartite hoppings and imaginary hoppings between 
the same sublattices and preserves the particle-hole sym- 
metry, this wavefunction (q;|$) can be written as a prod- 
uct of two Slater determinants AiDet{Mij)^ , where Ai 
is just an unimportant Marshall sign factor, and: 

= {[^A{h) + i'B{h)] 

+ (-!)* [i^Aih) - ^B{h)]}e"''-^^ (11) 

Here ijjA{k){')pB{k)) is the wavefunction amplitude on 
sublatticc A(B), rj is the coordinates of the up spins in 
configuration a, and ki is the momcntums in the momen- 
tum space. The Rcnyi entropy S2 of this wavefunction 
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can be calculated by VMC method detailed in the last 
section. 

Since TEE 7 is a subleading term in EE, naively it is 
rather difficult to extract since the leading boundary law 
term dominates the total entropy. However, as shown by 
Kitaev and PreskilU^, one can divide the total system 
into four regions A, B,C, D and consider the following 
linear combination of entropies that completely cancels 
out the leading non-universal term: 



S2,A + 82,3 
S2,ABC 



S2,C — S2,AB — S2,AC 



S2,BC + 

(12) 



where 82,11 is the Rcnyi entropy S2 corresponding to 
the region TZ. Note that the above formula neglects 

terms of 0(e~ «") where ^ is the correlation length and 
I' = min(lA,lB,lc)- This is because the above combi- 
nation of entropies cancels out all local contributions to 
the entanglement entropy which are expected to form 
a taylor series expansion in the variable following 
the arguments in Ref.— . One could still have terms 
that are non-perturbative in ^/l' which do not cancel 
out, hence the exponential dependence of the error made. 
This has also been confirmed by a direct calculation on 
select models in Refi^. Therefore, to accurately deter- 
mine TEE 7, it is important that ^ <C I' ■ Since we have 
a variational parameter A/t ar our disposal, one may 
tune it to minimze the correlation length ^. Chosing an 
optimal value A ~ 0.5t, and dividing the total system 
with dimensions 12 x 12 lattice spacings into La x La 
squares A and B, and an La x 2La rectangle C (see 
FigE}, one finds 7 = 0.343 ± 0.012 for an = 3 sys- 
tem and 7 = 0.344 ± 0.043 for an = 4 system^. 
Both are in excellent consistence with the expectation 
of 7 = log ("\/2) ~ 0.347 for its ground states' two fold 
degeneracy. 

To confirm the role of Gutzwiller projection in obtain- 
ing the QSL, one may also study the BCS wavefunction 
used to obtain CSL without Gutzwiller projection. This 
is an exactly solvable problem and one indeed finds a 
value w for TEE. This also serves as a benchmark for 
the VMC calculation. For an La = 3 system, the VMC 
calculation gives 7 = — 0.0008±0.0059, in agreement with 
vanishing TEE^. 

One may also study the effect of increasing correlation 
length on the TEE results. As mentioned above, as cor- 
relation length increases, he finite size effects from sub- 
leading terms become more important. Indeed, as seen in 
FigHl the TEE calculated using VMC deviates monoton- 
ically from the from the universal value of 7 = log {V^) 
as we lower the gap size A/t for a system with typical 
length scale La = 3. 

As mentioned above, the CSL ground state can be writ- 
ten ~ Det{AlijY where Mij is the Slater determinant 
corresponding to a Chern insulator. This leads to gen- 
eralization to wavefunctions that are cube or the fourth 
power of a Chern insulator wavefunction. For example, 
consider the wavefunction: 
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Figure 3: The separation of the system into subsystem A, B, 
C and environment, periodic boundary condition is employed 
in both X and y directions. 
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Figure 4: Illustration of finite size effect: chiral SL TEE 7 as 
a function of the relative gap size given by 2 A/t for character- 
istic system length La = 3. The larger the relative gap, the 
closer the ideal value is approached. For comparison, 7 for 
chiral SL at La = 4, and a large relative gap, is also shown. 
On the same plot, TEE of a lattice version of = 1/3 Laugh- 
lin state at La = 3, 2A/t = 1.0 is also shown. The dashed 
lines are the ideal TEE values of 7 = log (%/2) for the chiral 
SL and log (\/3) for the v = 1/^ Laughlin state. 



*i/3(''i, r2 . . . , rjv) = ^^(ri, . . . , r^) 



(13) 



where $ is the Chern insulator Slater determinant de- 
fined above. Clearly, the product is a fermionic wave- 
function, since exchanging a pair of particles leads to a 
change of sign. This is similar in spirit to constructing 
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the corresponding Laughlin liquid of m = 3 of fermions, 
by taking the cube of the Slater determinant wavefunc- 
tion in the lowest Landau level 

iP{zi,..., zn) ^'[[{z^- Zj)e '"^ (14) 

i<j 

However, unlike the canonical Laughlin state, composed 
of lowest Landau level states, these are lattice wave- 
functions. An interesting question is whether the low- 
est Landau level structure is important in constructing 
states with the topological order of the Laughlin state, or 
whether bands with identical Chern number is sufficient, 
as suggested by field theoretic arguments. 

Calculating the TEE using VMC with La ^ 3 one 
finds 7 = 0.5894 ± 0.0272 for the m ~ 3 wavefunction, 
consistent with the ideal value 7 = log (V^) = 0.549. 
Similarly, for the fourth power of the Chern insulator 
Slater determinant, one finds 7 = 0.732 ± 0.076, again 
consistent with ideal value in the thermodynamic limit: 
7 = log (^/4) = 0.693. 

These results offer direct support for the TEE formula 
7 = log D as well as their validity as topological ground 
state wavefunctions carrying fractional charge and statis- 
tics. The lattice fractional Quantum Hall wavefunctions 
discussed here may be relevant to the recent studies 
of flat band Hamiltonians with fractional quantum Hall 
states^--^. 



B. Beyond TEE: Extracting Quasiparticle 
Statistics Using Entanglement 

At the fundamental level, topologically ordered phases 
in two dimensions display a number of unique proper- 
ties such as nontrivial statistics of emergent excitations. 
For example, in Abelian topological phases, exchange 
of identical excitations or taking one excitation around 
another (braiding) leads to characteristic phase factors, 
which implies that these excitations are neither bosonic 
nor fermionic. A further remarkable generalization of 
statistics occurs in non- Abelian phases where excitations 
introduce a degeneracy. An important and interesting 
question is whether the ground state directly encodes 
this information, and if so how one may access it? It 
is well known that topologically ordered phases feature a 
ground state degeneracy that depends on the topology of 
the space on which they are defined. Also, as discussed 
in the previous sections, the ground state wavefunctions 
of such states contain a topological contribution to the 
quantum entanglement, the topological entanglement en- 
tropy (TEE). In this section we show that combined to- 
gether, these two ground state properties can be used to 
extract the generalized statistics associated with excita- 
tions in these states. 

The generalized statistics of quasiparticles is formally 
captured by the modular S and U matrices, in both 
Abelian and non-Abelian states^"—. The element Sa 



of the modular S matrix determines the mutual statis- 
tics of i'th quasiparticle with respect to the j'th quasi- 
particle while the element Uu of (diagonal) 14 matrix de- 
termines the self-statistics ('topological spin') of the I'th 
quasiparticle. Note, these provide a nearly complete de- 
scription of a topologically ordered phase - for instance, 
fusion rules that dictate the outcome of bringing together 
a pair of quasiparticles, are determined from the modu- 
lar S matrix, by the Verlinde formula^^. Previously, Wen 
proposed^ using the non-Abelian Berry phase to extract 
statistics of quasiparticles. However, the idea in Ref. 
requires one to have access to an infinite set of ground- 
states labeled by a continuous parameter, and is difficult 
to implement. Recently, Bais et alj^ also discussed ex- 
tracting S matrix in numerical simulations, by explicit 
braiding of excitations. In contrast, here we will just 
use the set of ground states on a torus, to determine the 
braiding and fusing of gapped excitations. 



It is sometimes stated without qualification, that TEE 
is a quantity solely determined by the total quantum di- 
mension D of the underlying topological theory. How- 
ever, this holds true only when the boundary of the re- 
gion A consists of topologically trivial closed loops. If the 
boundary of region A is non-contractible, for example if 
one divides the torus into a pair of cylinders, generically 
the entanglement entropy is different for different ground 
states (see Figure Indeed as shown in Ref. for a 
class of topological states, the TEE depends on the par- 
ticular linear combination of the ground states when the 
boundary of region A contains non-contractible loops. 
We will exploit this dependence to extract information 
about the topological phase beyond the total quantum 
dimension D. 



Let us briefly summarize the key ideas involved 
in extracting braiding statistics using quantum 
entanglemenli^^ and then discuss the details. We 
recall that the number of ground states N on a torus 
corresponds to the number of distinct quasiparticle 
types^i^i^^. Intuitively, different ground states are 
generated by inserting appropriate fluxes 'inside' the 
cycle of the torus, which is only detected by loops 
circling the torus. We would like to express these 
quasiparticle states as a linear combination of ground 
states. A critical insight is that this can be done using 
the topological entanglement entropy for a region A that 
wraps around the relevant cycle of the torus. To achieve 
this, we introduce the notion of minimum entropy states 
(MESs) 'E.j, j = 1 — N, namely the ground states 
with minimal entanglement entropy (or maximal TEE, 
since the TEE always reduces the entropy) for a given 
bipartition. These states are generated by insertion of a 
definite quasi-particle into the cycle enclosed by region 
A. With this in hand, one can readily access the modular 
S and U matrices as we now describe. 
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Figure 6: Vectors wi and W2 that define a lattice with peri- 
odic boundary conditions. That is, any point f on the lattice 
satisfies r = f + wi = f + W2. The area of the lattice is 
X TO2I where x denotes the cross product. 




Figure 5: Two types of entanglement bipartitions on the 
torus: (a) A trivial bipartition with contractible boundaries 
for which the TEE 7 = log D, and (b) A bipartition with non- 
contractible boundaries, where the TEE depends on ground 
state. 



ordered states where rotations can induce a change in the 
quasiparticle type. For example, in the Wen-Plaquette 
model2i, the rotations convert a charge-type excitation 
to a flux-type excitation. In such topological states, the 
procedure outlined above can give incorrect results, and 
we exclude such states from our discussion. 



1. Extracting Statistics from Topological Entanglement 
Entropy 

As mentioned above, for Abelian phases, the ij'th en- 
try of the S matrix corresponds to the phase the I'th 
quasi-particle acquires when it encircles the j'th quasi- 
particle. The U matrix is diagonal and the ii'th entry 
corresponds to the phase the i'th quasi-particle acquires 
when it is exchanged with an identical one. Since the 
MESs are the eigenstates of the nonlocal operators de- 
fined on the entanglement cut, the MESs are the canoni- 
cal basis for defining S and U. The modular matrices are 
just certain unitary transformations of the MES basis. 
For example. 



Sal: 



1 

D 



(15) 



Here D is the total quantum dimension and x and y 
are two directions on a torus. Eqn IlSI is just a unitary 
transformation between the particle states along different 
directions. In the case of a system with square geome- 
try, the S matrix acts as a 7r/2 rotation on the MES 

basis In general, however, x and y do not need 

to be geometrically orthogonal, and the system does not 
need to be rotationally symmetric, as long as the loops 

defining and inter wind with each other. We 

should mention one caveat with Eqn IlSI When we write 
Egn llSi we assume that the quasiparticlcs transform triv- 
ially under rotations. However, there exist topologically 



More concretely, let us denote the set of ground state 
wavefunctions as |^o,) where a = 1 — N and N is the 
total ground state degeneracy. To obtain the modular 5" 
matrix, one calculates the minimum entropy states |Si „) 
and |S2,q) for entanglement cuts along the directions wi 
and W2 respectively. Doing so leads one to the matrices 
Ui and U2 that relates the states |<^q) to |Si and |S2^q) 
respectively. The modular S matrix is given by 1/2^ Ui 
upto an undetermined phase for each MES. The existence 
of an identity particle that obtains trivial phase encircling 
any quasi-particle helps to fix the relative phase between 
different MESs, requiring the entries of the first row and 
column to be real and positive. This completely defines 
the modular S matrix. Therefore, the modular S matrix 
can be derived even without any presumed symmetry of 
the given wave functions. 

The above algorithm is able to extract the modular 
transformation matrix S and hence braiding and mu- 
tual statistics of quasi-particle excitations just using the 
ground-state wave functions as an input. Further, there 
is no loss of generality for non- Abelian phases, which can 
be dealt by enforcing the orthogonality condition in step 
2 which guarantees that one obtains states with quan- 
tum dimensions da in an increasing order. As shown in 
Ref. [tI, if the lattice has 27r/3 rotation symmetry, then 
the operation of rotating MES states by an angle of 27r/3 
corresponds to the action of US on MES. Therfore, in 
the presence of such a symmetry, one can extract the U 
matrix as well. 
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2. Exam-pie: Semionic Statistics in Chiral Spin-Liquid from 
Entanglement Entropy 

Let us revisit Chiral Spin-Liquid (CSL), first discussed 
in the Sec IV Al Even though topological properties of 
CSL are well established using field-theoretic methods^, 
from a wavefunction point of view, CSL cannot be dealt 
with analytically unlike continuum Laughlin states. For 
example, the low energy theory of CSL predicts a parti- 
cle with a semionic statistics (that is, the particle picks 
up a factor of -1 when it goes 27r around an identical 
one) which is difficult to verify in the lattice wavefunc- 
tion directly. Wc will show that the entanglement based 
algorithm discussed in the last subsection readily demon- 
strates the existence of semion as an excitation in CSL. 

As discussed in Sec IV A[ CSL corresponds to a d-|-id- 
wave state for chargeless spinous. The corresponding 
ground state wavefunction is obtained by Gutzwiller pro- 
jecting the mean-field BCS ground state wavefunction. 
CSL has a topological ground-state degeneracy of two 
on torus. One of the ground states can be obtained by 
enforcing periodic boundary condition along i-direction 
and anti-periodic along ij. We denote it by |0,7r). The 
orthogonal ground state is given by \tt, 0) with the same 
notation. 

Our task is to obtain the modular iS-matrix since it 
encodes the mutual statistics of quasiparticles. To do so, 
as outlined in the previous subsection, we consider linear 
combinations of |0,7r) and |7r, 0): 



cos 



0|O,7r) +sin(?!)|7r,0) 



(16) 



and numerically minimize the entanglement entropy with 
respect to the parameter so as to obtain the MESs^^. 
For system with square geometry the S matrix describes 
the action of 7r/2 rotation on the MESs. Since the two 
states |0,7r) and |7r,0) transform into each other under 
7r/2 rotation, this implies that S matrix is given by 
V^axV where V is the unitary matrix that rotates the 
{|0,7r) |7r,0)} basis to the MESs basis. From Fig El one 
finds that 
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Figure 7: The separation of the system into subsystem A, B, 
C and environment, periodic or antiperiodic boundary condi- 
tion is employed in both x and y directions, a: The subsys- 
tem ABC is an isolated square and the measured TEE has 
no ground state dependence, b: The subsystem ABC takes 
a non-trivial cylindrical geometry and wraps around the y 
direction, and TEE may possess ground stated dependence. 



Even though the S matrix obtained using our method is 
approximate, the quasiparticle statistics can be extracted 
exactly. In particular, the above S matrix tells us that 
the quasi-particle corresponding to dp = 1 docs not ac- 
quire any phase when it goes around any other particle 
and corresponds to an identity particle as expected, while 
the quasi-particle corresponding to di/2 = 1 has semion 
statistics since it acquires a phase of tt when it encir- 
cles another identical particle. Numerical improvements 
can further reduce the error in pinpointing the MES and 
thereby leading to a more accurate value of the S matrix. 

It is also interesting to compare the full dependence of 
TEE on the angle (j). Theoretically, one finds the follow- 
ing expression for TEE^^: 



27-7' 



log 



3 + sin (40) 



(18) 



where 7 is TEE for a region with contractiblc boundary 
while 7' is that for a region with non-contractible bound- 
ary. From FiglSl we sec that the numerically evaluated 
TEE is in rather good agreement with the exact expres- 
sion. 



/ cos(0.147r) -sin(0.147r)e"^ \ 
sin(0.147r) cos(0.147r)e''^ ) 

The existence of an identity particle requires positive 
real entries in the first row and column and implies = 0, 
which gives: 

f 0.77 0.63 \ _ J_f 1-08 0.89 \ 
0.63 -0.77 ) ^ ^\ 0.89 -1.08 ) 

which is in good agreement with the exact result 



5exact-^(^| -1) ^^^^ 



VI. GAPLESS QSLS AND ENTANGLEMENT 

In this section we consider Gutzwiller projected Fermi 
liquid wave functions which are considered good ansatz 
for ground states of critical spin-liquids^. We analyze 
two different classes of critical spin-liquids: states that 
at the slave-particle mean-field level have a full Fermi 
surface of spinous and those with only nodal spinous. 
The basic idea wc use here is that for free fermions, the 
boundary law is violated by a multiplicative logarithm, 
that is, S{1) ~ l'^-^ \og{l)^^. This suggests that for 
critical spin-liquids, whose low-energy theory contains 
a spinon Fermi surface^i, a similar result might be ex- 
pected. Indeed, one finds that for such Gutzwiller pro- 



10 




-0.2 -0.1 0.0 0.1 



Figure 8: The black dots show the numericaUy measured 
TEE 27 — 7' for a CSL ground state from hnear combina- 
tion 1$) = cos |0, tt) + sm4>\-K,Q) as a function of (j) with 
VMC simulation*^ using geometry in Fig. [TJa • Here 7 is 
TEE for a region with contractible boundary while 7' is that 
for a region with non-contractible boundary. The solid curve 
is the theoretical value from Eqn. (TH] The periodicity is ti/2. 
The red dots show the TEE for the same linear combination 
for a trivial bipartition. In the latter case, TEE is essentially 
independent of <j) and again agrees rather well with the theo- 
retical expectation (the dashed red curve). 



jected wavcfunctions the entanglement entropy scales as 
S{1) ~ l\og{l) for small system sizcs^. On the other 
hand, for wavefunctions with nodal spinons, the low- 
energy theory is believed to a two-dimensional CFT for 
which the result in Eqn|4]must hold. Thus, in this case, 
the VMC algorithm can be used to extract the universal 
shape dependent subleading constant. 

The wave-functions for these spin-liquids are con- 
structed by starting with a system of spin-1/2 
fermionic spinons fra hopping on a finite lattice at 
half-filling with a mean field Hamiltonian:i7jv/j;^ = 
'Ylirr' \.~^'rr' frafr'a + h.c?\^. The Spin wave-fuuctiou is 
given by 10) = Pg\4>)mf where \4>)mf is the ground 
state of Hmf and the Gutzwiller projector Pq = 
(1 — rii^ni^) ensures exactly one fermion per site. The 
one dimensional case was discussed in Ref. l76l 



A. Critical QSL with Spinon Fermi Surface 

Consider the mean-field ansatz that describes spinons 
hopping on a triangular lattice. In this case, one might 
expect that the projected wavefunction describes the 
ground state of a QSL whose low energy theory contains 
a spinon Fermi surfac e^ Indeed, one finds^ that for 
a triangular lattice of total size 18 x 18 on a torus and 
a subregion A of square geometry with linear size La up 
to 8 sites, the Renyi entropy scales as La fog^A (Fig. 
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Figure 9: Renyi entropy data for the projected and unpro- 
jected Fermi sea state on the triangular lattice of size 18 x 18 
with La = 1 — 8. Note, projection barely modifies the slope, 
pointing to a Fermi surface in the spin wavefunction. 



|9]). This is rather striking since the wave-function is a 
spin wave-function and therefore could also be written in 
terms of hard-core bosons. This result strongly suggests 
the presence of an underlying spinon Fermi surface. In 
fact the coefficient of the La logi^i term is rather similar 
before and after projection. This observation may be ra- 
tionalized as follows. Consider one dimensional fermions 
with flavors, where the central charge (and hence the 
coefficient of the log La entropy term) are proportional to 
N . However, Gutzwiller projection to single occupancy 
reduces the central charge to TV — 1. This is a negligible 
change for large A'^. Similarly, a two dimensional Fermi 
surface may be considered a collection of many indepen- 
dent one dimensional systems in momentum space, and 
Gutzwiller projection removes only one degree of free- 
dom. 

The area-law violation of the Renyi entropy for 
Gutzwiller projected wave-functions substantiates the 
theoretical expectation that an underlying Fermi surface 
is present in the spin wavefunction. 



B. Critical QSL with Nodal Spinons 

Next, consider a mean- field ansatz where the spinons 
hop on a square lattice with a flux tt through every 
plaquette^l. Thus, at the mean-field level, one obtains 
spinons with Dirac dispersion around two nodes, say, 
(7r/2,7r/2) and (7r/2,— 7r/2) (the location of the nodes 
depend on the gauge one uses to enforce the tt flux). 
The projected wave-function has been proposed in the 
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" Projected Triangular 




amounts to taking correlations into account, one might 
wonder whether the Fermi surface undergoes a mag- 
netic instabihty after the projection. Indeed, one finds 
that the non-zero magnetic order in the projected wave- 
function^, consistent with an independent study Ref. [tI. 
This was verified by calculating spin-spin correlations on 
a 42 X 42 lattice. Rcnyi entropy results for a a lattice 
of total size 24 x 24 with region A being a square up 
to size 12 X 12 are shown in the Fig. [TOl Though it is 
difficult to rule out presence of a partial Fermi surface 
from Renyi entropy, there is a significant reduction in 
the Renyi entropy as well as the coefficient of L^ZogL^ 
term as compared to the unprojected Fermi sea. 



Figure 10: Comparison of the Renyi entropy data for the pro- 
jected Fermi sea state on the triangular lattice and square 
(with and without 7r-flux) lattice as a function of the perime- 
ter P of the subsystem A. Here C is the constant part of 
the 82- We find S2 ~ Plog{P) + C for the projected trian- 
gular lattice state while S2 ~ P + C for the projected vr-flux 
square lattice state. For the square lattice state (no fiux), the 
projection leads to a significant reduction in S2 and the data 
suggests at most a very weak violation of the area-law in the 
projected state. 

past as the ground state of an algebraic spin liquid. The 
algebraic spin-liquid is believed to be describable by a 
strongly coupled conformal field theory of Dirac spinous 
coupled to a non-compact SU{2) gauge fieldi^iS. Be- 
cause of this the algebraic spin-liquid has algebraically 
decaying spin-spin correlations which can be verified ex- 
plicitly using Variational Monte Carlo^. 

Square lattice being bipartite, the projected wavcfunc- 
tion satisfies Marshall's sign rule and hence one is able to 
perform Monte Carlo calculation of Renyi entropy on big- 
ger lattice sizes compared to the triangular lattice case. 
Consider an overall geometry of a torus of size La x 4L^, 
La < 14, with both region A and its complement of sizes 
L^x2L^ (the total boundary size being La+La = "^La)- 
One finds^ that the projected wavefunction follows an 
area law akin to its unprojected counterpart and has the 
scaling 5*2 ~ O.SILa + 7 where 7 1.13 is a universal 
constant that depends only on the aspect ratio of the 
geometry^S. 



C. Effect of Projection on Nested Fermi Surface 

Finally, consider a mean-field ansatz as the fermions 
on a square lattice with nearest-neighbor hopping. The 
unprojected Fermi surface is nested. Since projection 



VII. SYNOPSIS AND RECENT 
DEVELOPMENTS 

In this article, we reviewed the perspective provided by 
the concept of entanglement entropy on QSLs. In par- 
ticular, we discussed the practical uses of the concept of 
TEE to detect topological order in realistic QSLs and a 
detailed characterization of topological order using en- 
tanglement as well. Returning to the questions posed in 
the introduction, let us discuss whether entanglement en- 
tropy is of any value in finding new QSLs? Remarkably, 
one of the most accurate methods to study strongly corre- 
lated systems, namely. Density Matrix Renormalization 
Group (DMRG)^, employs entanglement entropy as a 
resource to find ground states of particular Hamiltoni- 
ans. In brief, DMRG implements a real-space renormal- 
ization group that successively coarse-grains the system 
on the basis of the entanglement between a sub-region 
and the rest of the system. Though DMRG is most 
accurate for one-dimensional quantum systems, recent 
works are now beginning to explore the realm of two- 
dimensional strongly correlated systems. As an output, 
DMRG yields an approximate ground state of the sys- 
tem along with the entanglement entropy for a bipar- 
tition of the system into two halves. Having access to 
the entanglement entropy implies that one can use the 
results reviewed in the previous sections, thus allowing 
one to establish (or rule out) a QSL ground state. As 
an example, using DMRG Yan et aliS found that the 
ground state of kagome lattice Heisenberg Hamiltonian 
is a QSL. In a more recent study, Jiang et alii established 
the topological nature of this spin-liquid by calculating 
the TEE of the ground state for a region A with non- 
contractible boundary and found it to be close to log(2), 
a value that is apparently consistent with that of a Z2 
QSL. Similar results are found for other models which 
are known to have topologically ordered ground statesii. 
However, as discussed in Sec lVB ll TEE for a region with 
non-contractible boundary can be any number between 
and log(2), the maximum value of TEE (hence min- 
imum value for total entropy), being attained when the 
ground state corresponds to quasiparticle states. These 
results indicate that DMRG might be biased towards a 
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linear combination of ground states that minimize total 
entropy. As discussed in Sec lVBi entanglement can be 
used to determine quasiparticle statistics from a topo- 
logically degenerate set of ground states. Recently, using 
such an approach, Cincio and Vidal^ extracted the mod- 
ular S and U matrices for a microscopic lattice Hamil- 
tonian that has a Chiral Spin Liquid ground state. En- 
tanglement based measures have therefore emerged as a 
powerful conceptual and numerical tool in the context 
of quantum spin liquids. One may expect that in fu- 
ture, the confluence of ideas of many body physics and 



quantum information will throw up other equally fruitful 
directions. 
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